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I .  INTRODUCTION 


The  eroding  rod  model  for  deep  penetration1.2  is  attractive  because 
of  its  simplicity  and  its  ability  to  make  qualitative  predictions  that 
appear  to  be  useful  for  parametric  studies.  Nevertheless,  it  contains 
several  obvious  flaws.  Chief  among  these  is  the  use  of  the  modified 
Bernoulli  equation  and  of  the  oversimplified  rigid/perfectly  plastic 
material  model,  which  is  implied  by  that  equation.  A  critical  review  of 
the  model  as  it  presently  exists  is  presented  in  the  next  section. 

In  spite  of  its  shortcomings  the  eroding  rod  model  appears  to  be  a 
good  starting  place  for  an  experimental  investigation  of  penetration, 
which  in  turn  should  lead  to  a  more  complete  model.  In  the  third  section 
of  this  paper  the  theoretical  framework  for  an  experimental  program  is 
described.  The  theory  of  one  dimensional  wave  propagation  is  used  to 
show  how  data  from  instrumented  long  rods  and  targets  may  be  fitted 
together  to  give  a  coherent  picture  of  the  time  sequence  of  events  during 
penetration.  Data  for  one  impact  condition  is  then  compared  within  the 
theoretical  framework.  In  the  final  section  the  results  to  date  are 
discussed. 


II.  REVIEW  OF  THE  ERODING  ROD  MODEL 

The  eroding  rod  model1’2  for  deep  penetration  by  long  rods  is  based 
on  three  simple  ideas.  First,  a  long  rod  penetrator  may  be  considered 
to  be  rigid  and  undeforming,  except  where  it  is  in  contact  with  the 
target.  There  the  material  erodes  and  flows  at  a  constant  characteristic 
stress.  Second,  there  is  a  constant  characteristic  flow  stress  in  the 
target.  Third,  a  modified  Bernoulli  equation  holds,  which  includes  the 
characteristic  rod  and  target  stresses.  The  rod  stress,  of  course,  acts 
to  decelerate  the  intact  rear  end  of  the  rod.  The  situation  is  shown  in 
Fig.  1.  The  equations  expressing  these  ideas  are  as  follows2. 


Balance  of  Mass: 


i  =  p  -  u 


Balance  of  Momentum: 
Modified  Bernoulli: 


p  Z  u  =  -  E 
P 

Z  +  h  Pp(u-p) 


CD 


i  -2 

T  +  h  Pt  p 

stress  at  rod/target  interface 


1V.  P.  Alekseevskii,  "Penetration  of  a  Rod  into  a  Target  at  High  Veloc¬ 
ity,  "  Comb . ,  Expl. ,  and  Shook  Waves ,  2 ,  99-106  (1966). 

ZA.  Tate ,  "A  Theory  for  the  deceleration  of  Long  Rods  After  Impact ,"  J. 
Mech.  Rhys.  Solids ,  15,  387-399  (1967),  and  "Further  Results  in  the 
Theory  of  Long  Rod  Penetration,"  J.  Mech.  Phys.  Solids,  17,  141-150 
(1969). 
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Figure  1.  Long  Rod  Penetrating  a  Target. 


Here  £  is  the  instantaneous  length  of  intact  rod  and  u  is  its  instantan¬ 
eous  speed.  The  depth  of  penetration  is  p  and  the  rate  of  penetration 
is  p.  The  densities  of  rod  and  target  respectively  are  p  and  p^..  The 
characteristic  stresses  are  E  and  T  in  the  rod  and  target? 


Equation  (l)g  is  solved  for  p  and  situated  in  (1)^.  This  procedure 

gives  three  ordinary  differential  equations  for  £,  p,  and  u,  which  are 
to  be  solved  simultaneously.  The  integration  is  not  quite  straight  for¬ 
ward,  however,  because  there  are  simple  physical  constraints  on  p,  namely 


0  <  p  <  u. 


(2) 


The  lower  limit  is  required  since  penetration  can  only  open  a  hole  in 
the  target,  never  close  one.  The  upper  limit  is  required  since  it  is 
assumed  that  penetrator  and  target  are  always  in  contact,  but  since  the 
penetrator  cannot  increase  in  length,  £  £  0  in  (l)j. 

Now  consider  the  two  limiting  cases  in  turn  as  was  done  by  Tate2 . 
When  p  =  0,  there  is  no  further  penetration  and  no  further  flow  of  target 
material,  so  the  stress  in  the  target  must  be  less  than  T.  If  the 
penetrator  is  still  flowing,  then  the  target  stress  at  the  interface 
from  (1)^  is  given  by 

x  =  I  +  %>pu2  <  T  (3) 


Clearly  this  can  only  occur  if 

T  >  Z 


and 

u  S 


2(T-Z)/p 


(4) 
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Similarly,  if  p  =  u,  then  z  =  0,  there  is  no  further  flow  of  rod  mater¬ 
ial,  and,  therefore,  the  stress  in  the  rod  must  be  less  than  I.  If  the 
target  material  is  still  flowing,  from  (1)^  the  rod  stress  at  the  inter¬ 
face  is  given  by 

a  =  T  +  hPt  u2  S  Z  (5) 

This  can  only  happen  if 


T  <  Z 


and 


u  * 


E 


2(Z  -T)/p 


(6) 


It  must  be  concluded  that  the  material  model  implicit  in  the  use  of  the 
modified  Bernoulli  equation  is  the  highly  idealized  case  of  the  rigid/ 
perfectly  plastic  solid  for  both  penetrator  and  target  as  shown  in 
Fig.  2. 
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Figure  2.  Rigid/Perfectly  Plastic  Material. 
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Equations  (1)  must  be  replaced  by  the  following  set  of  equations, 
which  satisfy  the  restrictions  outlined  above. 


Z  =  p  -  u 


u 


P 


-  (T  +  %Ptu2)/p pZ,  if 


(6)  holds 


-  E/p  Z,  in  all  other 
P 

0,  if  (4)  holds 
u,  if  (6)  holds 


cases 


in  all  other  cases. 


(7) 


These  equations  are  simple  enough  to  be  solved  on  a  programmable  hand 
calculator.  Typical  integrations  of  (7)  giving  depth  of  penetration  as 
a  function  of  impact  speed  in  non-dimensional  form  are  shown  in  Fig.  3. 
The  results  appear  to  be  qualitatively  correct,  especially  the  character¬ 
istic  S-curve  for  the  case  T>E,  as  shown  by  the  experimental  results  of 
Stilp  and  Hohler3  and  Tate1*.  Equations  (7),  therefore,  have  consider¬ 
able  intuitive  appeal,  combining  as  they  do  simplicity  with  qualitative 
accuracy,  but  they  are  difficult  to  use  quantitatively  because  values 
for  the  characteristic  stresses  E  and  T  are  not  readily  available  a 
priori.  The  usual  procedure  has  been  to  choose  the  stresses  so  as  to 
fit  an  experimental  S-curve,  and  then  to  consider  them  as  material  con¬ 
stants.  Due  to  the  approximations  inherent  in  the  model,  this  procedure 
is  unlikely  to  be  satisfactory.  The  difficulties  are  threefold. 

The  first  problem  is  in  the  origin  of  the  modified  Bernoulli  equa¬ 
tion  itself.  In  the  theory  of  perfect,  incompressible  fluids,  the 
Bernoulli  equation  is  obtained  in  steady  flow  by  integration  along  a 
streamline.  This  gives 


P 


+  h  pv 


2 


constant 


(8) 


3F.  Hohler j  and  A.  J.  Stilp3  "Penetration  of  Steel  ccnd  High  Density  Rods 
in  Semi-Infinite  Steel  Targets 3"  Third  Int.  Symp.  Ballistics, 3  Karlsruhe 3 
Germany  (1977). 

**A.  Tate3  K.  E.  B.  Green3  P.  G.  Chamberlain  and  R.  G.  Baker 3  "Model  Scale 
Experiments  on  Long  Rod  Penetration3  "  Fourth  Int.  Symp.  Ballistics 3 
Monterey 3  CA  (1978). 
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Figure  3.  Penetration  as  a  Function  of  Impact  Speed. 


where  P  is  pressure.  Then  by  applying  the  result  at  a  stagnation  point 
and  on  the  stagnation  streamline  far  from  the  stagnation  point  in  both 
penetrator  and  target,  and  by  appealing  to  the  idea  of  continuity  of 
forces  at  the  stagnation  point,  the  Bernoulli  equation  as  applied  to  jet 
interaction  is  derived5. 


^  pp  v 


H  P. 


w 


(9) 


The  speeds  v  and  w  are  measured  with  respect  to  the  stagnation  point. 

Now  apply  the  same  ideas  to  the  steady  motion  of  a  rigid/perfectly 
plastic  solid.  With  the  z-axis  and  centerline  of  the  rod  coinciding  in 
Fig.  1,  the  z-component  of  the  momentum  equation  may  be  written 


t 

zz,z 


+  t 


zy.y 


+  t  =  p 

zx,x  H 


(u  +  U  U  +  U  +UU 

z,t  x  z,x  z,y  zz, 


On  the  centerline  u  =  u  =  0  by  symmetry,  and  for  steady  motion  u  *  0. 

^  y  z ,  t 

Shear  stresses  also  vanish  on  the  centerline,  but  not  their  derivatives, 
so  integration  between  points  a  and  b  yields 


t 

zz 


00  -  t 


zz 


(.a) 


dz 


(10) 


Equation  (10)  is  similar  to  (8)  with  the  important  difference  of  the 
remaining  integral.  If  equation  (10)  is  to  apply  to  the  rod,  point  b 
could  be  fixed  in  the  rod/target  interface,  and  point  a  could  be  located 
at  the  rigid  plastic  boundary.  The  exact  location  of  that  boundary  is 
unknown,  but  within  the  terms  of  the  material  model  the  stress  there 

must  be  the  constant  flow  stress,  t  =  -£,  and  the  particle  speed  as  it 

z  z 

enters  the  region  of  steady  flow  must  be  the  speed  of  the  rigid  rear  end 
relative  to  the  interface,  u-p.  Similar  considerations  must  apply  to 
the  target  with  the  z-component  of  stress  at  the  rigid/plastic  boundary 
being  t  =  -T  and  the  particle  speed  relative  to  the  interface  being  p. 

Since  stress  must  be  continuous  at  the  interface  and  since  the  stagnation 
velocities  are  zero,  the  modified  Bernoulli  equation  should  read 


%.  I.  Gurevich 3  Theory  of  Jets  in  Ideal  FtuidSj,  Academic  Press 3  New 
York  and  London  (1965). 
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CUD 


where  I 


rb 

I  (t  +  t  )  dz 

J  a  zy,y  zx.x' 


Thus,  even  within  the  limitations  of  the  assumptions  of  a  rigid/plastic 
material  and. steady  flow,  the  modified  Bernoulli  equation  cannot  be 
strictly  true. 

The  second  problem  concerns  the  validity  of  those  assumptions ; 
namely  steady  flow  and  a  rigid/plastic  material.  Taken  overall,  kinetic 
energy  penetration  is  clearly  an  unsteady  process,  especially  in  the 
initiation  and  final  stopping  or  breakout  stages,  but  there  may  be  parts 
of  the  process,  localized  in  both  space  and  time,  that  are  nearly  steady, 
as  for  example,  a  local  interaction  region  near  the  penetrator /target 
interface  during  the  intermediate  stages  of  penetration.  In  such  a 
domain  equation  (11)  should  be  valid,  but  even  then  there  is  no  clear 
way  to  make  an  accurate  estimate  of  the  characteristic  stresses.  The 
stress/strain  curve  of  most  real  materials  is  not  well  approximated  by 
the  rigid/plastic  assumption  because  of  work  hardening.  For  example,  a 
typical  curve  is  shown  in  Fig.  4.  Because  the  failure  processes  are 


—  2 


€ 


Figure  4.  Stress/Strain  Curve  with  Work  Hardening 
and  Rigid/Plastic  Approximation. 
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unknown  there  is  no  way  a  priori  to  choose  the  characteristic  or  average 
stress  E  either  on  the  basis  of  energy  equivalence  or  maximum  stress  or 
any  other  way.  In  the  target  the  situation  is  even  more  muddled.  For 
the  rod  it  has  been  tacitly  assumed  that  the  state  of  stress  is  nearly 
uniaxial,  but  in  the  target  the  stress  will  be  triaxial  with  a  substan¬ 
tial,  but  unknown  spherical  component.  The  matrix  of  components  on  the 
centerline  may  be  represented  by 


T1 

0 

0  ) 

1  P 

0 

0 

-S  0  0 

0 

T2 

0  (  =  - 

1  0 

P 

0 

+ 

Ohs  0 

0 

L 

0 

T3J 

1  0 

0 

P 

0  Ohs 

l 

The  first  matrix  on  the  right  contains  the  pressure  and  the  second  con¬ 
tains  the  deviatoric  components,  which  are  related  to  the  instantaneous 
flow  stress,  F,  in  the  usual  theory  of  plasticity6  by 

S  =  |  F  (13) 

The  target  resistance  stress  T  is  given  by 

T  =  -Tl=P  +  §  F  (14) 


Because  of  work  hardening  and  unknown  failure  processes,  the  target  flow 
stress  is  as  ill  defined  as  the  rod  stress,  but  in  addition  there  is  an 
unknown  hydrostatic  pressure,  which  may  vary  greatly  during  sequential 
stages  of  penetration.  A  final  consequence  of  real,  non-rigid  material 
behavior  is  that  there  will  be  considerable  plastic  deformation  at  some 
distance  from  the  quasi-steady  interaction  zone.  In  the  rod  this  has 
the  effect  of  slowing  material  down  (relative  to  the  rear  end)  before  it 
reaches  that  zone  so  that  the  speed  u  to  be  used  in  (11)  or  (1)  is  less 

than  the  speed  of  the  rear  end.  How  much  less  depends  on  the  rate  of 
work  hardening.  In  the  target,  plastic  deformation  has  the  effect  of 
increasing  the  spatial  rate  of  penetration:  that  is  to  say,  since  real 
target  materials  can  move  ahead  of  the  penetrator,  the  interface  between 
penetrator  and  target  will  move  farther  in  space  than  it  would  if  the 
target  were  rigid  in  non-plastic  regions.  Furthermore,  inertia  will 
tend  to  increase  the  crater  size  even  after  the  active  driving  forces 
have  ceased. 

The  third  problem  is  that  the  stress  tending  to  decelerate  the  rod 
in  (1)2  is  an  average  stress  over  the  whole  cross  section,  whereas  the 

stress  that  enters  into  the  modified  Bernoulli  equation  is  the  local 
stress  on  the  centerline.  These  may  be  somewhat  different  from  each 
other. 
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In  summary,  then,  it  has  been  shown  that: 

(1)  The  modified  Bernoulli  equation  cannot  hold  rigorously  because 
of  shear  stress  gradients  in  solids; 

(2)  The  deformation  processes  are  steady  at  most  in  a  limited 
domain,  but  even  then  it  is  not  possible  to  choose  the  characteristic 
stresses  £  and  T  a  priori  because  of  work  hardening  and  unknown  failure 
processes.  Furthermore,  plastic  deformation  will  extend  well  outside 
the  quasi-steady  zone;  and 

(3)  The  local  centerline  stress  is  unequal  to  the  average  cross- 
sectional  stress. 

For  these  reasons  the  eroding  rod  model  is  difficult  to  use  quanti¬ 
tatively  and  gives  only  limited  insight  into  the  details  of  the  actual 
rod/target  interaction,  but  in  spite  of  its  shortcomings  the  model  appears 
to  contain  the  germ  of  a  sound  theory  of  penetration,  as  shown  by  the 
qualitative  success  of  its  predictions.  The  model  probably  averages 
over  the  forces  and  erosion  processes  in  some  way  so  as  to  achieve  par¬ 
tial  success. 


III.  THEORETICAL  FRAMEWORK  FOR  AN  EXPERIMENTAL  PROGRAM 

An  experimental  program  has  been  devised  to  begin  probing  the 
details  of  the  interaction  process.  The  program  is  based  on  the  follow¬ 
ing  assumptions  and  theoretical  considerations.  Penetration  is  controlled 
by  flow  and  failure  processes  that  occur  at  or  near  the  penetrator/target 
interface.  These  processes  will  depend  on  the  material  properties  of 
the  particular  materials  involved,  and  therefore,  they  will  depend  only 
on  the  local  stress  and  strain  fields.  These  ideas  are  at  least  compati¬ 
ble  with  the  formation  of  a  region  of  steady  flow  at  the  interface. 

Since  it  is  a  local  phenomenon,  it  seems  reasonable  to  suppose  that  the 
steady  conditions  established  will  depend  only  on  the  local  geometry  and 
the  mass  flux  through  the  region.  In  particular  it  should  not  depend  at 
all  on  the  length  of  the  penetrator.  This  suggests  that  a  long  rod  with 
strain  gages  attached  could  be  used  as  a  probe  into  the  target  interior. 

If  the  rod  is  long  enough  so  that  wave  reflections  from  the  rear  cannot 
interfere  with  measurements,  then  the  interpretation  of  the  strain  data 
will  be  considerably  simplified.  At  the  same  time,  any  inferences  con¬ 
cerning  the  forces,  strains,  or  particle  velocities  at  the  impact  end 
will  be  equally  valid  for  shorter  rods  of  the  same  material.  In  other 
words,  it  seems  plausible  that  material  characteristics  will  fix  the 
boundary  conditions  locally  at  the  impact  end.  It  is  the  function  of 
the  experimental  program  to  determine  what  those  boundary  conditions  are. 

In  normal  impact  the  motion  of  the  rod  is  axi-symmetric.  If  the 
rod  is  treated  as  a  one  dimensional  continuum  where  radial  motions  are 
ignored,  the  equations  of  motion  may  be  written  as  follows: 
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(15) 


a 


x 


pu 


t 


In  (15)  a  is  the  axial  stress  in  the  rod  (referred  to  the  initial  cross- 
sectional  area),  p  is  the  initial  density,  u  is  the  particle  velocity, 
and  e  is  the  engineering  strain.  Subscripts  denote  partial  differentia¬ 
tion  where  x  is  the  material  coordinate  in  the  rod  and  t  is  time.  The 
spatial  location,  x,  of  a  material  particle,  X,  is  denoted 


x  =  x(X,t) 


we  have 


u  =  xt,  e  =  xx  -  1 


(16) 


(17) 


If  strain  gages  are  located  at  multiple  stations  along  the  rod  in 
an  impact  experiment,  then  strain  as  a  function  of  time  can  be  recorded 
at  each  station,  and  by  interpolation  strain  will  also  be  known  at  all 
other  stations  as  well.  Instrumented  rod  experiments  of  this  type  have 
been  conducted  by  Hauver  over  the  past  several  years.  He  uses  a  light- 
gas  gun  where  an  evacuated  target  chamber  and  accurate  alignment  permit 
data  of  unusually  high  quality  to  be  obtained.  Complete  experimental 
details  are  reported  by  him  in  other  papers7’8,  but  the  experimental 
conditions,  in  brief,  were  as  follows.  A  long  steel  rod,  (S-7  tool 
steel,  hardened  to  RC47) ,  with  length  of  254mm  and  diameter  of  8.13mm, 
impacts  at  normal  incidence  against  a  25.4mm  plate  of  rolled  homogeneous 
armor  at  RC27.  Speed  of  impact  was  nominally  1000  m/s  for  most  tests, 
and  710  m/s  in  a  few  others.  Strain  was  measured  in  the  rod  as  a  func¬ 
tion  of  time,  usually  at  stations  20,  40,  60,  and  80mm  from  the  impact 
end. 


With  e-t  data  at  hand,  curves  of  constant  strain  may  be  plotted  in 
X-t  space.  Along  such  a  curve  the  slope,  defined  as  Cp,  the  plastic 

wave  speed,  may  be  measured  and  related  to  derivatives  of  strain  as 
follows . 


7G.  E.  Hauver,  "Penetration  with  Instrumented  Rods, "  Int.  J.  Eng.  Sai. , 
16,  871-877  (1978). 

8£.  E.  Hauver,  "Experiments  with  Instrumented  Long-Rod  Penetrators, " 
Fifth  Int.  Symp.  Ballistics,  Toulouse,  France  (1980). 
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(18) 


de  =  e  dt  +  e  dX  =  0 

t  X 

dX  ft 

3T  -  P  ‘  "  «x 


These  curves  need  not  be  straight  lines,  but  if  they  are,  as  all  experi¬ 
ments  to  date  indicate,  the  analysis  is  greatly  simplified  because  cp 

then  depends  only  on  e.  Equations  (15)  may  be  integrated  with  respect 
to  X  along  lines  t  =  constant  as  was  described  by  Kolskyy  and  others. 

If  the  rod  has  speed  uq  and  zero  stress  and  strain  at  any  station  before 
a  wave  arrives ,  we  have 


u(X,t) 


a(X,t)  = 


fX 

uo  ♦  I  Ettx,t5 

J  00 

u  -  ( 

J  o 

I  e  Vx’ 


dX 


e(X,t) 

c  (e)  de 
P 


t)  dX 


(X,t)  2 

J  PCp(£)  de 


(19) 


(20) 


Note  that  the  stress  has  been  obtained  in  (20)  without  recourse  to  any 
constitutive  assumption.  In  fact,  a  dynamic  stress-strain  relation  has 
been  computed  from  the  experimental  data  using  only  the  simple  one 
dimensional  equations  of  motion  and  compatibility. 

Equations  (19)  and  (20)  determine  the  particle  velocity  and  stress 
everywhere  in  the  rod  in  principle.  Thus,  if  the  rate  of  rod  erosion 
were  known,  these  quantities  would  be  known  for  the  rod  as  it  enters  the 
quasi-steady  region  near  the  penetrator/target  interface.  The  situation 
is  shown  in  Fig.  5.  The  rays  of  constant  strain  are  taken  from  Ref.  8. 

The  timing  of  gage  failure  cannot  be  used  to  establish  the  unknown 
erosion  trajectory  because  gages  are  typically  destroyed  by  ejecta  before 
reaching  the  interface8 .  Independent  experiments  by  Netherwood  have 
measured  the  rate  of  penetration  into  the  target  directly  for  the  same 
materials  and  impact  conditions  as  in  Hauver's  experiments.  Complete 
details  are  given  in  his  paper,  but  essentially  his  technique  was  to 


SH.  Kolsky ,  Stress  Waves  in  Solids,  Dover  Publications ,  New  York  (1963). 

10P.  H.  Netherwood ,  "Rate  of  Penetration  Measurements,"  ARBRL-MR- 0297 8 > 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  Dec.  1979. 
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Figure  5.  Lines  of  Constant  Strain  in  X-t  Space  with 
Unknown  Erosion  Trajectory  [after  Hauver8]. 
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insert  contact  switches  into  the  target  block  through  holes  drilled  from 
the  side.  The  timing  of  switch  triggering  then  gives  the  trajectory  of 
the  penetrator/target  interface  in  target  material  coordinates.  Fig.  6 
shows  a  compilation  of  data  obtained  by  Netherwood10.  The  unsteady  entry 
region  is  clearly  shown,  followed  by  a  steady  intermediate  region  and  an 
erratic,  final,  unsteady  region. 

The  final  apparent  increase  in  speed  is  probably  due  to  pins  shear¬ 
ing  ahead  of  the  penetrator  as  a  small  plug  forms.  The  penetration 
trajectory  may  be  compared  with  strain  data,  measured  in  material  coor¬ 
dinate  in  the  rod,  if  both  sets  of  data  can  be  transformed  into  a  common 
coordinate  system. 

Fortunately,  all  data  may  be  easily  transformed,  at  least  approxi¬ 
mately,  into  laboratory  coordinates.  For  the  penetration  data  the 
transformation  is  a  simple,  time  dependent  translation.  As  penetration 
progresses  the  rear  surface  of  the  target  bulges  and  finally  fails  by 
plugging.  Netherwood  has  measured  this  motion  of  the  rear  surface  with 
a  streak  camera11.  If  deformation  between  switch  pin  and  rear  surface 
is  neglected,  then  the  translation  of  the  pin  in  front  of  the  penetrator 
is  the  same  as  that  of  the  rear  target  surface.  Target  deformation 
remains  completely  unknown,  but  would  contribute  only  a  small  correction 
compared  to  overall  translation. 

Rod  data  can  also  be  transformed  into  laboratory  coordinates.  Since 
particle  velocity  is  known  as  a  function  of  position  and  time  from  (19) , 
integration  with  respect  to  time  will  give  trajectories  in  laboratory 
coordinates  for  each  rod  station.  Hauver  has  computed  some  of  these 
trajectories8.  It  turns  out  to  be  more  instructive  to  plot  spatial 
trajectories  of  constant  strain.  From  (16),  with  X  =  X(t)  being  a  curve 
of  constant  strain  in  material  coordinates,  differentiation  yields  the 
following  equation. 


dx  3x  +  3x  dX 
dt  “  3t  3X  dt 

With  the  aid  of  (17)  ,  (18)  ,  and  (19)  this  becomes 

$£■  =  u  -  f  c  de  +  (1+  e) c  (21) 

dt  o  JQ  p  P 

Several  of  these  curves  are  plotted  in  Fig.  7.  Both  X  and  x  are  posi¬ 
tive  to  the  right,  and  uq  is  positive,  but  both  c^  and  e  are  negative 

since  strains  are  compressive,  and  plastic  waves  progress  to  the  left 
into  the  rod.  Note  that  curves  of  constant  strain  are  straight  lines  in 


11P.  H.  Netherwood ,  private  eommunieation. 
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Figure  6.  Envelopes  of  all  Data  Showing  Penetration  vs.  Time 
in  Target  Coordinates  from  Ref.  10. 


Figure  7.  Trajectories  in  Laboratory  Coordinates  of  Rod 
Stations,  Penetration,  Target  Free  Surface, 
and  Curves  of  Constant  Strain. 
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x-t  space  if  they  are  straight  in  X-t  space.  Also  shown  in  Fig.  7  is 
the  back  surface  motion  from  Ref.  11  and  the  translated  penetration 
trajectory.  In  some  of  his  experiments  Netherwood  recovered  the  residual 
penetrator,  whose  initial,  undeformed  length  was  determined  by  weighing, 
and  a  small  plug  approximately  5mm  thick.  The  trajectories  of  the  lead¬ 
ing  station  of  the  residual  penetrator,  labeled  96.3mm,  and  of  the  final 
plug  are  also  shown. 


IV.  DISCUSSION  OF  RESULTS 

Figures  6  and  7  contain  a  great  deal  of  useful,  albeit  incomplete, 
information  about  the  penetration  process.  In  Fig.  6  the  unsteady  entry 
region  seems  to  be  about  one  penetrator  diameter  in  width.  The  interme¬ 
diate  region  shows  a  nearly  constant  speed  of  penetration,  which  is  con¬ 
sistent  with  the  idea  of  a  zone  of  steady  state  deformation  in  penetra¬ 
tor  and  target. 

There  is  no  way  of  determining  whether  or  not  a  switch  triggers 
slightly  before  penetrator  arrival  in  material  coordinates,  but  since  it 
could  not  be  late,  the  measured  arrival  times  are  lower  bounds.  In  any 
case,  in  the  steady  central  section  of  penetration,  it  seems  reasonable 
to  assume  that  each  switch  is  early  by  the  same  amount,  on  average,  so 
that  the  true  curve  would  be  parallel  to  the  one  drawn  in  its  middle 
section  with  corresponding  minor  adjustments  in  the  unsteady  end  sections. 
The  width  of  the  breakout  zone,  about  9mm,  is  also  approximately  one 
penetrator  diameter,  although  it  is  suspicious  that  one  switch  pin  was 
usually  located  at  the  beginning  of  that  zone  and  may  have  played  a  role 
in  its  initiation.  Since  the  recovered  plug  was  only  about  5mm  thick 
rather  than  9mm,  it  seems  clear  that  erosion  of  penetrator  and  target 
continues  during  plug  formation  and  acceleration. 

In  Fig.  7  curves  of  constant  strain,  the  penetration  curve,  back 
surface  motion,  and  trajectories  of  selected  rod  stations  are  all  shown 
in  laboratory  coordinates.  There  is  actually  considerable  extrapolation 
of  data  here  since  strain  gage  records  end  near  the  target  surface  due 
to  interference  from  ejecta8,  and  the  penetration  curve  ends  at  the  start 
of  plug  formation.  The  free  surface  data  do  extend  to  150ys,  however, 
and  serve  to  anchor  the  trajectory  extrapolations  for  penetration  and 
for  the  96.3mm  station.  These  three  curves  (i.e.,  free  surface,  penetra¬ 
tion,  and  96.3mm  station)  together  with  the  known  plug  thickness  give  a 
remarkably  coherent  picture. 

Strains  higher  than  14%  were  measured  at  one  or  two  stations,  but 
determination  of  c^  becomes  unreliable  so  that  trajectories  for  higher 

strains  are  not  shown.  Nevertheless,  the  figure  clearly  suggests  that 
the  leading  edge  of  the  fan  of  constant  strains  might  be  roughly  parallel 
to  the  penetration  trajectory.  If  this  were  true,  then  equation  (21) 
would  give  an  estimate  of  rate  of  penetration  possible  for  various  levels 
of  maximum  strain.  At  the  very  least  it  will  give  a  lower  bound  estimate 
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for  p.  This  bound  as  a  function  of  strain  is  plotted  in  Fig.  8  for  the 
case  considered.  The  measured  value  for  p,  taken  from  Fig.  7,  is  about 
190  m/s,  which  is  not  much  greater  than  the  maximum  of  the  lower  bound. 
Note  that  the  bound,  as  given  by  equation  (21),  is  linear  in  the  impact 
speed  u  so  that  other  cases  may  be  considered  simply  by  moving  the 
horizontal  axis  up  or  down,  provided  the  shape  of  the  curve  does  not 
change  too  much  with  u  . 


Successful  measurements  have  now  been  made  of  I,  u,  and  p  in  equa¬ 
tion  (11) .  Another  series  of  experiments  to  measure  T  directly  for  the 
same  impact  conditions  are  also  under  way  by  Pritchard  at  BRL.  This 
leaves  only  the  two  terms  I  and  1^,  which  cannot  be  measured  directly, 

but  perhaps  could  be  estimated  analytically. 


In  future  experimental  work,  the  speed  of  impact,  target  thickness, 
and  penetrator  material  will  be  varied.  At  the  same  time  the  analysis 
is  being  extended  to  include  the  effects  of  radial  motion  of  rod  mater¬ 
ial  since  the  assumption  of  one-dimensional  motion  in  equation  (15)  is 
probably  not  completely  adequate. 


The  investigation  of  detailed  interactions  during  penetration  has 
been  guided  by  a  theoretical  framework  that  is  based  on  the  one 
dimensional  theory  of  wave  propagation.  Although  this  work  has  not  yet 
led  to  a  replacement  for  the  eroding  rod  model,  it  has  led  to  an  experi¬ 
mental  description  of  events  against  which  any  successor  model  may  be 
evaluated. 
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PARTICLE  SPEED 
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Figure  8-  Particle  Speed,  u,  and  Speed  of  Strain 
Wave,  dx/dt,  vs-  Strain,  e. 
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